Planar Spin Network Coherent States 



I. General Properties 

Donald E. Neville 
Department of Physics 

Temple University 
Philadelphia 19122, Pa. 

July 6, 2008 

Abstract 

This paper is the first of a series of three which construct coherent 
states for spin networks with planar symmetry. Constructing such 
states is not straightforward. After gauge-fixing, the full SU(2) sym- 
metry is broken to U(l); but one cannot simply use the U(l) limit 
of SU(2) coherent states, because the planar states exhibit an un- 
expected 0(3) symmetry arising from the closed loop character of 
the transverse directions. This paper uses an intuitive, rather than 
mathematically rigorous approach to construct a candidate set of 
0(3) symmetric coherent states. However, paper 2 of the series 
then demonstrates explicitly that the proposed coherent states are 
approximate eigenvectors of the holonomy and momentum operators 
E} (as expected for coherent states), up to small correction terms. 
These coherent states are superpositions of holonomies which obey 
the residual U(l) gauge symmetry only on average; that is, some 
holonomies in the superposition violate the symmetry, although the 
U(l) quantum numbers of these holonomies are peaked at values 
which obey the symmetry. An appendix discusses the closely re- 
lated case of cylindrical symmetry. 
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I Introduction 



This paper is based on the canonical, spin network formulation of 
quantum gravity. This is an approach built on a connection, rather 
than metric formulation of gravity [T[ The connection is real 
(3]. Basic variables are a densitized inverse triad and an exponen- 
tiated connection, i. e. a holonomy. Other operators of the theory 
are constructed by commuting holonomy with volume or Euclidean 
Hamiltonian The space on which the holonomy exists is a one- 
dimensional network of edges meeting at vertices, rather than a 
continuous manifold ^ E] . 

At present there is no empirical confirmation for the existence 
of this microscopic spin network structure, although in the future 
it may be possible to detect Planck-scale modifications to elemen- 
tary particle decay rates, or modifications to light propagation over 
cosmological distances [3 [3 [9], [10] . For now, the approach must be 
checked by undertaking calculations which confirm the internal con- 
sistency of the formulation, or confirm consistency with established 
principles. For example, area and volume operators for the theory 
possess a discrete spectrum |Jjy, in qualitative agreement with the 
conclusions from thought experiments that lengths smaller than a 
Planck length are not measurable [12]. Cosmological calculations, 
extended back to the big bang, yield a finite result [13] . Black hole 
calculations predict area oc entropy and yield a formula for the en- 
tropy [gl[T5]. 

One topic which requires further investigation is the classical limit 
of the spin network approach. The systems studied so far (black 
holes, homogeneous cosmologies) have reasonable classical limits, 
but they are so highly symmetric that they cannot propagate gravi- 
tational waves. In order for the system to propagate waves, the spin 
network Hamiltonian must be nonlocal: a single term in the Hamil- 
tonian must be able to change the spin network at two or more 
neighboring vertices [16]. Once the Hamiltonian is made nonlocal, 
it is no longer obvious that the constraints have the correct commu- 
tation relations, in the limit of fields varying slowly over many spin 
network vertices [17]. I.e., it is not obvious that the theory possesses 
general coordinate invariance in the classical limit. 

Several approaches use the path integral, or "spin foam" ap- 
proach to study the classical limit, rather than the canonical, spin 
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network approach. Aleschi and Rovelli calculate the gravitational 
Green's function, then check for a correct long-range behavior [T8] . 
This approach puts limits on allowable spin foam vertices. Amb- 
jorn, Jurkiewicz, and Loll put the system in a heat bath and solve 
numerically for the ground state, to see if the geometry is reason- 
able [in] . This work puts limits on the topology of admissible paths, 
and favors spin networks which are causal. Markopoulou, and also 
Oeckl, have studied the application of the renormalization group to 
gravity theorv[20ll2T]. 

Thiemann and Winkler develop another approach to the classi- 
cal limit: within the canonical approach, construct coherent states 
which can be used to study the theory approximately, in the classical 
limit [221E31I21]- This is the approach followed in the present paper, 
which constructs coherent states for a space possessing two commut- 
ing spacelike Killing vectors which may be written d/dx and d/dy. 
This planar symmetry is the simplest which allows propagation of 
gravitational waves, and therefore requires a nonlocal Hamiltonian. 

For this case one cannot use the coherent states constructed by 
Thiemann and Winkler. In one sense these states are too general. 
They possess the full SU(2) local gauge invariance, whereas planar 
symmetry allows the gauge to be fixed, until only U(l) gauge rota- 
tions around the Z axis survive [25] . (Lower case coordinates x,y,z,t 
refer to coordinates on the global manifold; upper case coordinates 
X,Y,Z,T refer to coordinates in local Lorentz frames.) In another 
sense, the Thiemann Winkler states are not quite general enough, 
because the planar states must possess an 0(3) symmetry which is 
not a limit of the SU(2) gauge symmetry [261 127] . 

This is the first of a series of three papers which construct coher- 
ent states for the planar case [281129]. In any discussion of coherent 
states, there is always a danger that one can get buried in detail. 
The discussion can degenerate into a collection of statements, with 
each one-sentence statement followed by a two-page proof. In an 
effort to avoid this, I have adopted a three-paper approach which 
allows readers to take in as much or as little detail as suits their 
purposes, as follows. 

The present paper gives a general overview and relatively few 
details. I will refer to this paper as paper 1. Section [III of the 
present paper describes a suitable Hilbert space for plane waves, 
and discusses the need for the 0(3) symmetry. This section is a 
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summary of an earlier paper by the author [26], and is included to 
make the present paper reasonably self-contained. 

Section IIIII constructs coherent states for the planar case. Thie- 
mann and Winkler derive their SU(2) coherent states in two ways. 
They use results of Hall on the mathematics of coherent states on 
group manifolds [3D]; and they show that Hall's methods parallel a 
procedure in standard quantum mechanics for turning a free parti- 
cle wavefunction into a minimal uncertainty wave packet. The latter 
derivation is perhaps easier to grasp intuitively, and I have followed 
it in section IIIII 

I follow an intuitive approach when constructing coherent states 
in the present paper. Therefore I must show explicitly that the 
states have the properties expected of coherent states. This is done 
in paper 2, the second paper in the series [2S]- This paper shows that 
the coherent states are approximate eigenf unctions of the operators 
which will be used to construct the Hamiltonian: 

O I coherent state) =< O >| coherent state) + SC, 

where the operator O is either a densitized inverse triad Ef or a spin 
one-half holonomy. The main body of paper 2 derives the leading 
contribution to the matrix element (the c- number function < O >) 
and gives qualitative estimates for the size of the small correction 
(SC) terms. Quantitative estimates of the SC terms are worked out 
in the appendices to paper 2. 

Some readers may not wish to contemplate the full detail of paper 
2, or they may prefer a short introduction before reading the details. 
For those readers, section HVl of the present paper gives an overview 
of the main steps in the calculation of matrix elements. The focus 
is on the structure of the calculation (which steps are done in what 
order, and why). This section contains a little detail, but not yet 
full detail. 

Cylindrical symmetry is very closely related to planar symmetry. 
Both metrics admit two commuting, spacelike Killing vectors. Ap- 
pendix [X] discusses the modifications needed to handle cylindrical 
rather than planar symmetry. 

Constructing a spin network Hamiltonian is a little like attempt- 
ing to navigate through a darkened room full of furniture. In the 
case of the room, it is virtually impossible to proceed in any direc- 
tion without barking one's shins against an obstacle, which requires 
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a choice of new direction. In the case of the Hamiltonian, one can- 
not proceed far without bumping into some ambiguity, which again 
requires a choice. 

The present series of papers might appear to be only an exphca- 
tion of the properties of coherent states, independent of any choice 
of Hamiltonian; however, one chooses to emphasize certain proper- 
ties of coherent states, because these are the properties one expects 
will be relevant later when calculating the Hamiltonian. There is at 
least an implicit choice of Hamiltonian involved in a paper of this 
type. I will try to make explicit relevant ambiguities and choices. 

One such choice is the assumption, made above, that the Hamil- 
tonian is constructed from spin 1/2 holonomies in the fundamen- 
tal representation of SU(2). This seems to contradict results from 
the spin foam approach to quantum gravity. In this approach, one 
starts with an initial spin network living in n space dimensions, then 
propagates it forward in time by gluing appropriate n+1 dimensional 
simplices to the initial network. For a concrete, relatively simple ex- 
ample, consider the Euclidean version of 2-1-1 gravity, with all edges 
labeled by SU(2) spins. The spin network version of this theory is 
given by Rovelli [31j- The edges which are added at each time step 
can have arbitrary spin; there is no limitation to new edges with 
spin 1/2 or to new spins which differ from old byj ± 1/2. The 
evolution does not seem to be generated by a Hamiltonian which 
singles out spin 1/2. 

Ideally, this particular shin-bark would dictate a unique change 
of direction: given the spin foam unitary evolution operator for 
2-1-1 gravity, discover the Hamiltonian which generates the evolu- 
tion. Presumably, the Hamiltonian would contain holonomies in 
representations of higher spin than 1/2. The natural choice of uni- 
tary evolution operator for this theory is a sum over a set of basic 
spin rearrangements (the so-called Pachner moves); but the Hamil- 
tonian which generates this operator has not been found [32]. (If 
the Hamiltonian itself is made into a sum of Pachner moves, it does 
not behave like a constraint; it does not annihilate physical states.) 
Spin foam theories which are more sophisticated than 2-1-1 gravity 
can be even harder to reconcile with the canonical approach [33] . 

At the present time, there seems to be a disconnect between the 
spin foam and canonical approaches. The spin foam approach can 
construct a convincing unitary evolution operator, but no Hamilto- 
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nian constraint. The canonical approach has a constraint, but the 
formahsm is frozen in time. In this situation, it is necessary to make 
a choice. I remain on the canonical side of the canonical-covariant 
divide and assume the holonomies are spin 1/2. By moving in this 
direction, of course, I may end up merely knocking over a lot more 
furniture. 

Although spin 1/2 is the simplest choice, I should report that 
Gaul and Rovelli have investigated Hamiltonians involving higher 
spin [31]. They find there is no problem in principle with using the 
higher holonomies. 

A third paper in the series contains some speculations about the 
form of the Hamiltonian [29]. Although paper 3 does use some re- 
sults from papers 1 and 2, it relies primarily on dimensional analysis 
arguments. In logical terms it is independent of the first two papers 
and should be evaluated independently. 

The coherent states constructed in this paper are not exact eigen- 
states of the residual U(l). Exact eigenstates could be obtained by 
angle-averaging the coherent states which will be constructed in sec- 
tion mil But these coherent states are already quite complicated, 
even before an angle average. It seems simplest to test for general 
covariance and gravitational wave propagation, initially using states 
which obey U(l) only on average; later, one can refine the calcula- 
tion. For studies of angle- averaged states which obey SU(2) exactly, 
see reference [55] . 

II The Planar Hilbert Space 

In reference [26] I constructed a kinematic basis for the planar case, 
and for completeness I review the highlights of that construction 
here. Call the direction of propagation the z direction. I will use 
lower case x, y, z, a, b,- ■ ■ for global coordinates, and upper case X, 
Y, Z, A, B,- ■ ■ for local SU(2) indices. 

Gauge fixing may be used to simplify the E and connection fields 
|25] . The off-diagonal elements E^ and E\, with a = x,y and A 
= X,Y, can be gauged to zero; similarly, A^ and A\ may be set 
to zero. This means that the holonomies along the longitudinal z 
direction are quite simple, involving only the rotation generator Sz 
for rotations around Z, while the transverse holonomies (those along 
the X and y directions of the spin network) involve Sx, Sy and are 
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rotations in the X,Y plane. Each transverse holonomy therefore has 
an axis of rotation of the general form 

m = (cos 0, sin 0,0), (1) 

for some angle 0. (More precisely, there is one holonomy for each 
transverse direction and one for each transverse direction, (f)^ and 
(py. Since the two directions are treated equally, I will discuss only 
the X holonomies, and will suppress the subscript x for now.) The 
corresponding spin 1/2 holonomy is 

= exp[-icTz(0 - 7r/2)/2] exp[icTy^/2] exp[+i(Tz(0 - 7r/2)/2] 
= /i(i/^)(-(/. + 7r/2,^,0-7r/2) (2) 

On the last two lines I have written the usual Euler angle decompo- 
sition for this rotation. 

Given the above, the natural basis for the transverse holonomies 
might seem to be the set of rotation matrices 

h(j)(-0 + vr/2,^,0-7r/2), 

where j is the highest weight obtained by multiplying together 2j 
/i(i/2) matrices. However, this basis is not convenient because it 
has complicated behavior under the action of the E . To obtain 
objects with simpler behavior, note that the matrix /i^^/^^ has only 
three independent components. ( Because the matrix is special, a 
rotation around an XY axis, the two diagonal elements are equal: 
= h ) Group the three independent components as follows: 

{Vx^Vy.Vz) 

= (sin(^/2) cos(0 - 7r/2), sin(^/2) sin(0 - 7r/2), cos(^/2)); (3) 
the inverse relation is 

h(i/2) = iVg + iV- (a- X Z), (4) 
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where boldface denotes a 2x2 matrix. Under action of an E , the 
three V components rotate infinitesimally hke the components of a 
three- vector. 

_ E^exp[i j A^Ssdx] 

E^VW = Vr,{N\SA\M){jK/2). (5) 

Subscripts A,B = X, Y only. On the second line there is a generator 
Sa on both sides of the holonomy, because the holonomy in a trans- 
verse direction is a closed loop, beginning and ending at the same 
vertex; therefore the E (which acts at the vertex) grasps the holon- 
omy at both ends. This grasp-at-both-ends feature, characteristic 
of closed loops, gives rise to the 0(3) symmetry. That symmetry 
may be exhibited by introducing the linear combinations V. From 
the last line, above, the grasp is just an infinitesimal 0(3) rotation, 
when expressed in terms of the Vs. 

Some technical points: In eq. I have omitted the delta func- 
tions coming from the [E ,A] commutators, because the delta func- 
tions are always canceled by the area and line integrals associated 
with E^ and A^. (I always suppress the area integration associ- 
ated with each E .) The overall factor 1/2 is a relic of the integrals 
over the deltas, 1/2 because the deltas occur at the endpoints of 
the edge integration, is the product of Immirzi parameter times 
SttG; h = c=l. 

I now return to the main point: because the V's transform more 
simply than the h's under the action of E , one should multiply 
together V's rather than h's to get a basis set of holonomies. Since 
the grasps generate an 0(3), and the V's form a vector under the 
action of this 0(3), the products of V's generate the usual spherical 
harmonics of 0(3). In particular, components of V itself form the 
simplest spherical harmonic of 0(3), 

{V±,Vo) = (Tsin(^/2)exp[±z(0-7r/2)]/v^,cos(e/2)); 
VM{h) = v/V(2LTl)lL=i,A/(^/2,0-7r/2). (6) 

Polynomials of degree L in components of V can be expanded in 
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a series of the higher harmonics Ylm, which therefore form a ba- 
sis. (More precisely, there are two bases, Ilx.Mx and Yiy^My for 
holonomies along the x and y directions respectively.) These har- 
monics transform simply under the action of E : 

(7/€/2)-1E^Ylm = EnYln(L,N I S± | L,M), (7) 

where Y^m = Yim{0/2, — 7r/2). The unconventional half-angle re- 
minds us of the origin of these objects in a holonomy /i^/^ depending 
on half-angles. 

The Y's are known to be proportional to matrix elements of ro- 
tations, 

v/47r/(2L + 1)Ylm = ^fl{-<P + vr/2, 0/2, - 7r/2). (8) 

Therefore eq. ([7]) is also correct if D's are substituted for Y's. I prefer 
D's to Y's in what follows, because use of D's (will require awkward 
factors of \/ ^ti j {2L + 1) in initial formulas, but) ultimately will 
result in fewer factors ^/4:Tc/ {2L + 1). The basic vector V{h)M is 
also simply related to a D: 

Vih)M = V47r/(2L + 1)Yl=i,m{0/2, (j) - 7r/2) 

= Dl^l{-^ + 71/2,9/2,^ -n/2). (9) 

From now on I will use a hat notation to emphasize that V is a unit 
vector. 

The planar calculation involves three groups: SU(2), because the 
holonomies in the Hamiltonian are rotation matrices in SU(2); 0(3), 
because the action of the grasps generates an 0(3) group; and U(l), 
because the usual SU(2) gauge invariance is broken to U(l) by the 
gauge fixing. It is worth taking a minute to contemplate when to 
use which group. 

Presumably the Hamiltonian will be constructed using holonomies 
/i(i/2) jj^ ^}^g fundamental representation of SU(2). The Euler de- 
composition of h^^^'^\ eq. ([2]), shows /i^^/^-* as depending on the full 
angle 9, but of course the actual matrix elements of h^^^'^^ contain 
half-angles sm{9/2),cos{9/2). When these matrix elements are re- 
arranged to form a vector D^^\h) in 0(3), the matrix elements of 
D^^\h) (and D^^\h), eq. ([8])) inherit this dependence on half-angles. 
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I am therefore using the same notation h for an SU(2) matrix 
and an 0(3) matrix; yet the SU(2) and 0(3) matrices have different 
dependence on the angle 0. Hopefully this will cause no confusion, 
because of the following circumstance. If the matrix has a super- 
script (1/2), as then it is in SU(2), and the Euler decompo- 
sition involves the full angle; if the matrix has a superscript other 
than (1/2), then it is in 0(3), and the Euler decomposition involves 
a half angle. 

As mentioned above, presumably the Hamiltonian will contain 
holonomies in the fundamental representation of SU(2), h^^/'^\ The 
coherent states, on the other hand, will be sums of representations 
of 0(3). Therefore when the Hamiltonian acts on the state, for 
calculational purposes one should replace the SU(2) /i^^/^^ matrices 
by the linear combinations D^^\h) which transform simply under 
0(3): 

+ 7r/2, e,(f)- 7t/2) I coh) 
^ Di^i,{-<P + 77/2, e/2, <p - 77/2) I coh) (10) 
= < Z^iif (-0 + ^/2, 0/2, - 77/2) >| coh) + SC, 

where the brackets denote the peak, or average value, and SC de- 
notes small corrections (typically down by 1/\/L). The details of 
the SU(2) — i> 0(3) replacement are given by eq. ([3]). Once the peak 
values of the D{l){h) are known, one can solve for the peak values 
of the /z^^/^) using the inverse equations, eq. (jl]). 

Now consider the relation between 0(3) and the residual gauge 
invariance U(l). Most of the 0(3) rotations have nothing to do 
with SU(2) gauge invariance. The exceptions are rotations around 
Z, which are identical to the residual U(l) gauge invariance. Proof: 
because of the gauge fixing, the axis of rotation for the matrices 
remain in the XY plane. This reduces the gauge rotations 
from SU(2) to the set of U(l) rotations around the Z axis. Under 
an infinitesimal Z rotation 6 is unchanged, while 0— >0 -|- 6(j). The 
corresponding changes in the matrices h^^^'^^ and the 0(3) harmonics 
are 
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The last line is identical to the result of an 0(3) rotation around 
the Z axis. □ 

It is straightforward to work out the consequences of U(l) gauge 
invariance. Each vertex has one x holonomy, one y holonomy, and 
two z holonomies (one entering or initial, one leaving or final). The 
total change in phase of the vertex under a Z rotation is therefore 

exp[i5/?(Mx + My + Mzf - Mzi)], 
where the Mz come from the Z holonomies, which go as 




U(l) invariance demands that the quantity 

(M^ + My + Mzf - Mzi) 

vanish, for each vertex. 

The coherent states constructed here will not have unique values 
for Mj. and My. The states are superpositions of D^^^hi) matrices 
(i = x,y); and the superposition contains a range of values Mj. The 
superposition is sharply peaked at central values of the Mj, however, 
so that all values of + My in the packet are very close to the 
average value 

< M, + My >=< S^"^ + Sl'^ > . (11) 

The above equation suggests we may have to study the matrix 
elements of spin, as well as E and holonomies, in order to study 
U(l) invariance. Fortunately, the E operators are identical to the 
spin operators, since E^ brings down a factor of S"^"*; we get spin 
matrix elements with no extra effort. 

Ill Coherent States 

This section proposes a set of candidate coherent states. I cannot 
use the coherent states available in the literature, since these are for 
the general case, full SU(2) symmetry, and do not exhibit the 0(3) 
symmetry discussed in the previous section. However, the results 
for full SU(2) motivate a tentative choice of states for the planar 
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case, together with a choice of parameters for labehng these states. 
The second paper in this series verifies in detail that the candidate 
states have the properties required of coherent states. 

What properties should a coherent state possess? The Hamilto- 
nian is constructed entirely from spin-half holonomies and E oper- 
ators. When these quantum operators act upon the coherent state, 
the result should be c-numbers. From the previous section, the 
holonomies should be thought of as spin one spherical harmonics 

dS(-0 + vr/2, ei2, - 7r/2) := D(i)(h)oM, 

while the E essentially multiply the state by a spin operator. A set 
of coherent states should therefore have the properties 

L,M 

D«(/^)oA I ...) = <Z^«(/^)oA>| + 
(7/t/2)-iE^ I . . .) = Y.DoNmL.N\SA\L,M)c{L,M-...) 

= <El >]...) + SC, (12) 

SC denotes small corrections, down by order l/y/< L >, < L > the 
peak value of L. Eq. (fT2|) means that the coherent state is an ap- 
proximate eigenvector of the holonomy and triad, with eigenvalues 
given by the numbers < D^^\h)QA > and < E^ >. If D(h) depends 
on angles {6, 0), then < D{h) > depends on angles {a, j3) which may 
be interpreted as the peak values of (6', 0). The dots (. . .) stand for 
coherent state labels to be introduced shortly. 

There are two kinds of SU(2) coherent states available in the 
literature. The first type has no sum over L in eq. fll2p . only a sum 
over M These coherent states are too simple for present 

purposes. If L is fixed, the E operators can be made classical, but 
not the holonomies: one must superimpose many L's to get sharply 
peaked angular values for (6*, 0). 

The second kind of coherent state was suggested by Hall |30] 
and elaborated for quantum loop gravity by Thiemann and Winkler 
[231 [21|- Their results were derived for the general case, full 
local SU(2) symmetry, whereas the planar case possesses only U(l) 
gauge symmetry. Simply taking a U(l) limit of the general case does 
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not work, since the planar limit must possess the 0(3) symmetry 
discussed in the previous section. 

However, the Thiemann Winkler coherent states may be under- 
stood intuitively as generalizations of the minimal uncertainty states 
for the free particle. One applies a certain recipe to construct the 
free particle states, and the same recipe works in the general SU(2) 
case. Once this intuitive approach is understood, it is straightfor- 
ward to use the same recipe to generate a set of candidate coherent 
states for the planar case. 

I now review the recipe for constructing a coherent state for the 
free particle. Start from a wave function which is a delta function. 



This wave function is certainly strongly peaked, but it is not nor- 
malizable. Also, it is peaked in position, but it should be peaked 
in both momentum and position. To make the packet normalizable, 
insert a Gaussian operator exp[— p^/ (2(T^)]. (Choosing the Gaussian 
form is a "cheat", because we know the answer; but for future ref- 
erence note that all the eigenvalues k"^ of must be positive, so 
that the Gaussian damps for all k.) To produce a peak in momen- 
tum, complexify the peak position: xq ^ xq -\- ipq/u"^. With these 
changes, the packet becomes 



The last hne follows after completing the square on the exponential, 
and exhibits the characteristic coherent state form. 

Note there is not just one state, but a family of coherent states, 
characterized by the parameter a. The shape of the wave function is 
highly sensitive to a; but the peak values Xq^Pq are independent of a, 
as is the minimal uncertainty relation AxAp — h/2. The coherent 
states constructed below contain a parameter t which is analogous 
to a. 




N / exp[-p7(2o-^)] exp[iA;(x - xq) + kpo/a'^] dk/2'K 
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Now apply the above recipe to the planar case. The position x 
becomes the pair of angular variables 0) on the group manifold. 
The complete set of plane waves become a complete set of spherical 
harmonics. To construct a delta function in angles, 

5{e/2 - «/2)(5(0 - sin(a/2) 

I introduce spherical harmonics D^^\u) depending on angles {a,P) 
in the same way that the D^^''{h) depend on {9, (p). 

D(^)(n)oM = D'^'^\-(3 + 7c/2,a/2,f3-n/2)oM 

= V47r/(2L + l)YLMia/2, P - n/2). (14) 

Compare eq. f|T^ to eq. ([8]). I can now write the delta function in 
angle as a sum over spherical harmonics. 

S{9/2-a/2)5{^-P)/sm{a/2) = 5^((2L+l)/47r)D(L)(h)oMD(L)(u)SM. 

L,M 

(15) 

As discussed in the last section, it is more convenient to use D 
matrices rather than Y^^m 's, but the reader who wishes to exhibit the 
latter can use eqs. (jSj) and (fill) to replace D's by Y's in eq. (ITB]) . The 
sum then takes on a form which may be more familiar, just ^ YY*. 
The wave number k in the free particle example corresponds to L in 
the planar case. 

Continue with the recipe: dampen the sum using a Gaussian 
exp[— tL(L + l)/2]. Complexify by extending the angles in u to 
complex values, replacing u by a matrix g in the complex extension 
of 0(3). The coherent state has the general form 

\u,p) = iV^((2L + l)/47r)exp[-tL(L + l)/2]D(^)(/i)oM/^^''H^7)oM 

L,M 

= NY,---D^''\h)oMD^'^\g])Mo. (16) 

L,M 

A vector p is needed to characterize the matrix g, as follows. 
Every matrix in SL(2,C), the complex extension of SU(2), can be 
decomposed into a product of a Hermitean matrix times a unitary 
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matrix ("polar decomposition"; see for example [38j)- E- g- for the 
fundamental representation, 

g = Hu = exp[cr ■ p/2]u. 

It follows that every matrix in 0(3) also has a polar decomposi- 
tion, obtained by restricting the representations of SU(2) to repre- 
sentations with integer spin. 

:= H^^^u^^^ (17) 

At this point I must choose six input parameters: the three Euler 
angles which determine the unitary matrix u and the three compo- 
nents of p, which define the complex extension. By analogy with the 
free particle case, if u determines the peak angles, then the complex 
extension p determines peak values of the canonically conjugate 
variable, the angular momentum. 

In this paper I use the following choices for these parameters. 
Restrict p to be an arbitrary vector in the XY plane. 

P3 = 0. (18) 

Restrict the unitary matrix u to be an arbitrary rotation with 
axis of rotation in the XY plane. The Euler decomposition of m^^^ 
is 

u(^) = u{-p + 7i/2,a,p -7i/2) 
:= exp[m ■ Sa/2]; 
n = (cos/3,sin/5,0). (19) 

(Compare this definition of u to the corresponding definition of the 
matrix h, given at eqs. ([1]) and ([2]). Both u and h have their axis of 
rotation in the XY plane.) 

These choices are certainly plausible. In the limit of no complex- 
ity (in the limit p — > 0) the coherent state will reduce to the original 
delta function and u will become the peak value of h. Since h has 
its axis of rotation in the XY plane, the peak value should also be 
a matrix with axis in the XY plane. 
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— * 

Also, the generator 5* ■ p for H is the complexification of the 
generator S ■ h for u. Since n hes in the XY plane, the same should 
be true for p. In short, since h is a very special matrix (its axis 
lies in the XY plane) the peak matrix u and its complexification H 
should also be special. 

In paper 2, I calculate the expectation values < > and 
< D'^^\h)QA > only for the case p^ = 0, and axis of u in the XY 
plane, eqs. ffTSl) and ffH?]) . At the end of section III of paper 2 I state 
a theorem which describes what happens when u and p are replaced 
by more general choices, 

(u,p)^(u,p), 

where p and the axis of u need not be in the XY plane. To keep an 
already too lengthy paper 2 as short as possible, I have not supplied 
a proof of the theorem; but the industrious reader can easily do so 
by using the techniques developed up to that point. Use of the more 
general values seems to add nothing but complexity. 

Since p now lies in the XY plane, p may be parameterized using a 
magnitude, p, plus an angle fi. It is convenient for later calculations 
to take /i to be the angle between p and n, the axis of u. 

p = p[cos(/? + /i),sin(/? + /i),0]. (20) 

(From eq. f|T9|) . (3 is the angle between n and the X axis.) 

When the expectation values < > and < D^^\h)QA > are 
calculated in paper 2, they turn out to be perpendicular. 

= ^<D«(h)oA><E^>. (21) 

A 

There is no constraint like this for the expectation values in the case 
of general SU(2) symmetry. 

The following quantum mechanical analogy is helpful (probably, 
essential) in clarifying the situation. Consider an electron moving 
in a central potential and described by polar and azimuthal coordi- 
nates {6/2, (/) — 7r/2). (The definitions of polar and azimuthal angles 
are unorthodox but acceptable.) The angular wavefunction for the 
electron, as well as the coherent state eq. f|T6|) . are both superposi- 
tions of spherical functions Y^m, or equivalently rotation matrices 
D{h)oM- Therefore eq. ffTB]) can serve as a coherent state describing 
the angular motion of an electron. 



16 



When the three quantities D^^\h)QA are written out, they turn 
out to be the three components of the (unit) radius vector of the 
electron. Therefore < D^^\h)oA > is the peak value of the unit ra- 
dius vector. Continuing with the analogy, < > is the peak value 
of the electron angular momentum vector, since the operator, 
when acting on a D(h), brings down a factor of spin Sa_- Because 
angular momentum is perpendicular to the orbit plane, these two 
vectors (< D^^\h)oA > and < E^ >) must be perpendicular. 

This constraint follows from the symmetry, not the dynamics (it 
does not matter whether the electron is moving in a hydrogen atom 
or in a spherical well); therefore the spin network coherent state 
must obey the same constraint. (The expectation values continue 
to obey the constraint even when the parameters u and p are gen- 
eralized to parameters {u,p) not lying in the XY plane.) 

IV Form of the Action on a Coherent State 

Paper 2 calculates the action of E^ and D^^\h)QA (equivalently, 
the holonomy /i^/^) on the coherent state | u,p). This section is 
intended to be an overview of those calculations. It summarizes the 
main steps in the two calculations, with a minimum of detail. 

Let Oa stand for either or D^^\h)oA', both have a vector index 
A. Then from eqs. (fT6|l and (fT7|) . 

Oa\u,p) = Ar^((2L + l)/47r)exp[-tL(L + l)/2] 

L,M 

■[OAD^'^\h)oM]D^''\u]H)MO 

= N J2 ((2L + l)/47r)exp[-tL(L + l)/2] 

L,M,L',M' 

■D'^^'\h)oM'{L'M' I Oa \ LM) 

■D^'^\u\)mnD^''\H)no- (22) 

(Neither E nor h changes the "zero" index on D^^\h)QM-) 

Eq. fl22|) is very general, but hard to interpret. However, since 
the actions of E^ and D(^)(/i)oa on the coherent state yield vectors, 
it will help to introduce an orthonormal triad of geometrically sig- 
nificant unit vectors and expand the right hand side of eq. fl22p in 
these vectors. 
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I use the rows of the orthogonal matrix D^^\u) as the orthonor- 
mal triad, since these provide a natural set of geometrically signifi- 
cant vectors. Proof: when Oa = D^^\h)oA in eq. ( l22l) (when Oa is 
a holonomy), one expects 

D^'\h)oA\u,p) = <D^'\h)oA>\u,p) + SC 
= D^'\u)oa\u,p) + SC; 

D^'\u)oA = Vih = u)A. 

(The last line is eq. ([9]).) I. e. the zeroth row of the matrix D^^\u) 
should be the peak value of the holonomy; and at least this row 
should occur when I expand the right hand side of eq. fl22l) . 

The remaining two rows, D^-^\u)±^a, are also geometrically sig- 
nificant. D^^\u)qa = V{u) is the image of the Z axis under the 
rotation u. Since u is a rotation around an axis n, it follows that 
h is perpendicular to V{u). The rows D^^\u)±^a are also perpen- 
dicular to V{u) = D^^\u)oA, therefore are linear combinations of 
n and the other unit vector perpendicular to V{u), namely n x V. 
It follows that the three rows of the orthogonal matrix D^^\u) are 
linear functions of the geometrically significant vectors n, V, and 
nxVn 

The exact relation between these three vectors and the rows of 
D^^\u) is worked out in paper 2. For now, note it is not too hard 
to produce a factor of D^^\u) in eq. (!22|) . Since the matrix ele- 
ment < Oa > in eq. (l22l) is essentially a Clebsch-Gordan coefficient, 
one can use the rotation properties of this coefficient to move the 
D^^\u]) matrix through the matrix element and produce factors of 
D«(m). 

{L'M' I Oa I LM)D^''\u])mn = D^'''\u])m'M"{L'M" \ Ob \ LN)D^'\u)ba; 
Oa\u,p) = N J2 ((2L + l)/47r)exp[-tL(L + l)/2] 

L,N,L',M' 

■D^'^'\hu\),M'{L'M' I Ob \ LN)D^^\u)baD^''\H)no- (23) 

The D'^'^\u)ba can be expressed in terms of n, V, and n x V", to 
reveal the component along each of these directions. 

Next, focus on the sum over L, N, L', M' in eq. fl2^ : how can 
it be simplified? We wish to prove that the coherent state | u, p) 
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is an (approximate) eigenstate of the operator Oa- This means it 
is desirable to simphfy the sum in eq. (123|) . until it resembles as 
much as possible the corresponding sum for the original coherent 
state. Eq. contains the product D^^'\hu])oM'D^^\H)m- The 
original coherent state, eq. (ITBi) . contains a similar product, but with 
indices M' and N replaced by a single index M, and L' = L. 

The indices M' and N in eq. (1231) are not necessarily equal. They 
differ by the index B, because of the selection rules obeyed by the 
Clebsch-Gordan coefficient in {L'M' \ Ob \ LN). 

N = M' -B. 

Also, when Oa is an E operator, L' will turn out to be L; but for 
Oa a holonomy, L and L' will differ by one unit: L' = L ± 1. 

Despite these mismatches, it is possible to rearrange the sum in 
eq. (123|) so that it closely resembles the sum in the original coherent 
state. The philosophy is to leave the D{hu]) factor alone; this factor 
is oscillatory, and hard to approximate. Work with the D(H) factor 
instead, which turns out to be exponential and relatively simple. 
Rewrite eq. fl23|) as 

Oa\u,p) = N J2 ((2L + l)/47r)exp[-tL(L + l)/2] 

L,L',M' 

■D^'^'\hu\)oM'D^'^\H)M'o 

■[{L'M' \Ob\L,M'- B)D^'^\H)M'-Bfll D^''\H)M'fl] 
■D^'\u)ba. (24) 
For Oa an E operator, L' = L, and the product 

D(L')(hut)oM'D(L)(H)M'o 

is now identical to the corresponding product in the original coherent 
state. For Oa a holonomy operator, L' = L ± 1; but it is best not 
to replace the D^^^{H)m'o factor by D^^'\H)m'q- The L' = L ±1 
states have peak values of L shifted by ±1 unit. Since this is a 
real physical effect which cannot be removed by changing the D(H) 
factors, it is simpler and clearer to leave the second line, eq. fl2^ . 
as is. 

We are not out of the woods yet: the square bracket in eq. ([21]) 
is a function of L and M, and cannot be taken out of the sum. 
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However, the factor D^^\H) mo ^^v[~tL(L + l)/2] turns out to have 
the normahzed Gaussian form typical of coherent states, 

D(^)(//)Moexp[-tL(L + l)/2] = 

>/tJ^exp[-t{L + 1/2 -p/t)V2] 

y/l/[{L + l/2)7r] exp{-MV[2(L + l/2)]}/(p, t). 

Therefore, to get the dominant contribution, one can expand the 
square bracket in eq. around the mean values, < L + 1/2 > = 
p/t and < M > = 0, and keep the first, constant term in the series. 
Let b(L, M; L'-L, B) denote the square bracket in eq. 1^^. 

b{L,M) = b{<L>,0) + {<L>dLb)AL/ <L> 
+(< L > dMb)AM/ < L > + ■ ■ ■ ; 
AL := L- < L>= L + 1/2 -p/t; 
AM := M- < M >= M. (25) 

I have suppressed the arguments L'-L and B, since they are held 
fixed. The expansion coefficients b(< L >, 0) and < L > dxb 
all turn out to have the same order of magnitude. Because of the 
Gaussian nature of D(H), the factors of AX/ < L > are order 
ax/ < L >, where ax is the standard deviation associated with the 
variable X (X = L or M). The first term, b(< L >, 0), is therefore 
the approximate eigenvalue of the operator Oa, while the remaining 
terms in the series are small corrections (SC) down by factors of 
ax/ < L>. 

The final results for the action of holonomy and triad on the 
coherent state are (with replaced by D^^l{h)) 

D;^}{h)\u,p) = [dI^I{u) = Va]\u,p) + SC; 
(7/t/2)""'^E^ \ u,p) = < L > {fiA cos fi — {n X V)a sin /i) \ u,p) + SC. 

The expectation value of spin Z component (relevant for U(l) con- 
servation; see eq. ffTTj) ) is 

< S^'^ + S^^ >=< Lx > sin(ax/2) sin/ix+ < Ly > sin(ay/2) sin/iy. 

This is the calculation in broad outline, and paper 2 fills in the 
details. 
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A note on the appearance of L+1/2 rather than L in Gaussian 
expansions of the D(H): the factor L+1/2 is an approximation for 

\/L(lTT) = L + 1/2 + order(l/L). 

This paper considers terms to order 1/\/L but does not always take 
into account terms down by 1/L. The factor of 1/2 is therefore an 
educated guess rather than the result of a calculation. The theory 
contains a number of scalars (the magnitude of S; standard devia- 
tions); and the scalar in SU(2) theory is \/ L{L + 1), rather than L. 
I use the easy-to-type L+1/2, rather than L, whenever it is likely 
that a scalar quantity should be \/L{L + 1), even though my ap- 
proximations may not be accurate enough to establish this. 

V Conclusion 

Note two nice features of the calculation sketched in section [IVl the 
use of the rotation properties of the Clebsch-Gordan coefficients at 
eq. (p3|) produces a convenient set of unit vectors; the subsequent 
power series expansion of the square bracket makes it easy to sepa- 
rate leading terms from small corrections, and evaluate the latter if 
needed. Both these features carry over to other calculations where 
the coherent states exhibit a continuous symmetry, including the 
general SU(2) calculation. 

The general SU(2) case involves a complete set of functions on a 
group manifold, whereas the present, 0(3) case involves a complete 
set of harmonics on a sphere. This seems to be a small difference; 
nevertheless, the 0(3) case is unusual in two respects. As shown 
in section IIIH the expectation values of E and holonomy must be 
perpendicular. Also, if peak values are tuned by means of an input 
matrix u and vector p, in the SU(2) case the matrix u always equals 
the peak value of holonomy h; whereas in the 0(3) case we shall 
see in paper II that u equals the peak value only for the specific 
parameter value ps = 0. 

I am not clear as to the deeper mathematical reasons for these 
differences. However, it seems clear that one should proceed with 
caution when trying to construct coherent states using a basis other 
than representations on a group manifold. 
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A Cylindrical Symmetry 



What happens when the symmetry is cyhndrical, rather than pla- 
nar? On first inspection, the cyhndrical case appears to be a simple 
relabeling of the planar case. The wave propagates in the r direc- 
tion in the cylindrical case, and a right-handed coordinate system 
is desirable. Therefore one can relabel 

X, y, z^0, z, r; X, Y, Z->$, Z, R 

The gauge fixing goes through as for the planar case [25]. In partic- 
ular, the E matrix is again block diagonal. The 1x1 block now con- 
tains E^, rather than E|; the 2x2 block now has rows and columns 
labeled by (0, z) and Z). Similarly, the matrix is also block 
diagonal. One can choose a representation of the generators such 
that "S'ij" is always diagonal (always given by Sz), while ''Si^, " and 
"S'z " are given by Sx and Sy respectively. Except for the relabel- 
ing, the cylindrical theory now looks exactly like the planar theory, 
suggesting that a solution to the former will also work for the latter. 

Can this be right? We know that in field theory the two symme- 
tries have different solutions, because the boundary conditions are 
different (different falloff with r or z, at r or z — >■ infinity; finiteness 
requirement at the origin, r = 0, for the cylindrical case). However, 
in a spin network formulation manifold indices, (r, 0, z) or (x, y, 
z), disappear. They are integrated over areas or curves in order to 
form spatial diffeomorphism invariants; and cylindrical and planar 
coordinates differ precisely by a spatial diffeomorphism. Somehow, 
in our haste to get rid of the coordinates on the manifold (ugly, 
because arbitrary), we have lost the distinction between cylindrical 
and planar. We still have the local Lorentz indices (R, $, Z) or (X, 
Y, Z); but these are too local to reveal the difference between the 
two symmetries. 

Perhaps the solution is to recover some global properties in a spin 
network context? Yes, but hiding the manifold coordinates was a big 
step forward, and one would like to avoid a step backward. It should 
be possible to introduce global properties without going back to the 
underlying manifold coordinates. E. g., for the cylindrical case, one 
must impose finiteness of the total energy within a unit length of 
cylinder; and for the planar case, finite total energy per unit area. 
Since the energy is given by the integral of a surface term at spatial 
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infinity, in field theory these finiteness requirements imply different 
asymptotic behaviors for cylindrical and planar fields. 

In spin network theory, one can require a finite total energy, with- 
out going to the underlying field theory or reintroducing manifold 
coordinates. The spin network is only superficially one-dimensional. 
Since both transverse holonomies are integrated over an unspecified 
repeat distance, the spin network surface term will be two dimen- 
sional, an integral over the two transverse repeat distances. It is 
not clear how large this transverse area is, but one can divide the 
energy integral by a transverse area integral. 



The resultant energy per unit area is independent of the unknown 
repeat distances. 

(I have imagined the transverse repeat distances as arbitrary, 
but fixed during the calculation. Alternatively, one can assume the 
symmetric theory is the limit as the repeat distance is taken to zero 
|39j . I am not sure which point of view is correct, or whether it 
matters. However, both points of view would seem to require the 
area integral in the denominator. 

In field theory, the linearly polarized cylindrical case is especially 
interesting because it can be solved exactly PU]. In the spin net- 
work approach, however, it may not be straightforward to impose 
the constraints which reduce the wave to a single linear polariza- 
tion ^27j. The field theoretic constraints set an E and its conjugate 
momentum equal to zero. In the spin network approach, the con- 
jugate momentum, a component of the extrinsic curvature, is con- 
structed by commuting the Euclidean Hamiltonian with the volume 
operator. Since the spin network Hamiltonian is non-local (since 
it involves holonomic loops connecting two or more vertices), the 
conjugate momentum becomes non-local. 

Furthermore, in order to solve the field theoretic case, one must 
(specialize to linear polarization, then) canonically transform to a 
new coordinate system: R, T coordinates, in the notation of Kuchaf 
|10] . The new T coordinate is a function of a momentum which 
Kuchaf calls vr^. In the notation of the present paper, that momen- 
tum is 




E| dx dy or / E^ d0 dz. 




S'K^E^. 
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K is the extrinsic curvature. The prime on the summation means 
omit r and R in the sum over indices. Again, this operator is not 
defined in a spin network context (or in quantum field theory for 
that matter); again, it must be defined using commutators of volume 
with Hamiltonian; again, it is non-local. 

This is as far as one can go without a specific Hamiltonian. In 
the field theoretic case, the linear polarization constraints may be 
used to simplify the Hamiltonian. Despite the non-locality, A similar 
simplification may occur in the spin network case. 
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